arXiv: 1506.0453Ivl [astro-ph.SR] 15Jun2015 


Mon. Not. R. Astron. Soc. 000, [T]-?? (2015) Printed 16 June 2015 (MN MgC style file v2.2) 


SCORPIO: 

A deep survey of Radio Emission from the stellar life-cycle 

G. Umana^*, C., Trigilio\ T.M.O. Franzen^, R.P. Norris,^, P. Leto^ A. Ingallinera^ 
C.S. Buemi^ C. Agliozzo^’^, F. Cavallaro^’^’^, L. Cerrigone^ 

^INAF- Osservatorio Astrofisico di Catania, Via S. Sofia 78, 95123, Catania, Italy 
^CSIRO Astronomy and Space Science, PO Box 76, Epping, NSW 1710, Australia 
^ Universitd di Catania, Dipartimento di Fisica e Astronomia, Via Santa Sofia, 64, 95123 Catania, Italy 
^Millennium Institute of Astrophysics, Santiago, Chile 
^ Universidad Andres Bello, Avda. Republica 252, Santiago, Chile 

^ASTRON, the Netherlands Institute for Radioastronomy, PO Box 2, 7990 TiA Dwingeloo, The Netherlands 


ABSTRACT 

Radio emission has been detected in a broad variety of stellar objects from all stages of 
stellar evolution. However, most of our knowledge originates from targeted observations of 
small samples, which are strongly biased to sources which are peculiar at other wavelengths. 
In order to tackle this problem we have conducted a deep 1.4 GHz survey by using the Aus¬ 
tralian Telescope Compact Array (ATCA), following the same observing setup as that used 
for the Australia Telescope Large Area Survey (ATLAS) project, this time choosing a region 
more appropriate for stellar work. In this paper, the SCORPIO project is presented as well as 
results from the pilot experiment. The achieved rms is about 30 pJy and the angular resolu¬ 
tion ^10 arcsec. About six hundred of point-like sources have been extracted just from the 
pilot field. A very small percentage of them are classified in SIMBAD or the NASA/IPAC 
Extragalactic Database (NED). About 80% of the extracted sources are reported in one of 
the inspected catalogues and 50% of them appears to belong to a reddened stellar/galactic 
population. The evaluation of extragalactic contaminants is very difficult without further in¬ 
vestigations. Interesting results have been obtained for extended radio sources that fall in the 
SCORPIO field. Many bubble-like structures have been found, some of which are classified at 
other wavelengths. However, for all of these sources, our project has provided us with images 
of unprecedented sensitivity and angular resolution at 2.1 GHz. 

Key words: Galaxy: stellar content - Radio continuum: stars - radio continuum: ISM - stars: 
evolution - stars: formation - techniques: interferometric. 


1 INTRODUCTION 

In the last few years, large radio surveys such as the NR AO VLA 
Sky Survey (NVSS) (ICondon et al.ll 19981) . the Faint Images of the 
Radio Sky a t Twenty-centimeters (FIRST) (iHelfand etal.|[i999h 
and ATLAS (iNorris et al.l l2006l : iHales et afl l20l4) have revealed 
different populations of radio emitting objects. However, very few 
stars have been found and nearly all known radio stars have been 
detected by targeted observations directed at small samples of stars 
thought likely to be radio emitters. The radio detection of differ¬ 
ent st ar types, covering different stages of stellar evolution (iGudell 
l2002h . strongly suffers from limited sensitivity and from selection 
bias as they have been obtained by targeted observations aimed to 
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address a specific problem related to a particular aspect of radio 
emission. 

This approach has proven to be quite productive but it is biased 
against discovering unknown, unexpected, or intrinsically rare ob¬ 
jects, preventing a good knowledge of radio stars at the sub-mJy 
level. In consequence, at the moment, it is quite difficult to provide 
a trustworthy forecast on the real potential of the Square Kilometer 
Array (SKA) and its pathfinders in the field of stellar radio astron¬ 
omy. From an anal ysis aimed to point o ut the stellar population in 
the FIRST survey, iHelfand et al.l (Il999l) concluded that even at a 
threshold of ^ 1 mJy stellar radio emission is quite rare. Similar 
results were found bv iKimball et al.l (l2009l) . These results are not 
surprising as FIRST and other radio deep surveys were designed 
for extragalactic studies and only regions at high Galactic latitude 
were covered. The space density of stars, unlike that of extragalac¬ 
tic sources, varies with Galactic latitude, being much more concen- 
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trated in the disk. To overcome this problem, it is clear that a deep 
blind survey, carried out in a sky patch well suited for stellar work, 
is necessary. 

1.1 Stellar radio emission 

The improvement of the observational capabilities have led to the 
discovery of radio emission in a broad variety of stellar objects 
from all stages of stellar evolution. In many cases, radio obser¬ 
vations have revealed astrophysical phenomena and stellar activ¬ 
ity not detectable by other means. Broadly speaking, the brightest 
stellar radio emission appears to be associated with active stellar 
phenomena such as flares, related to t he presence of a strong and/o r 
variable magnetic held, or mass-loss dSeaauistll 1 993l : lGudelll2002h . 
Much of our knowledge of microwave emission from radio stars 
comes from the study of active stars and binary systems as a large 
fraction of these have been found to be strong radio sources. Their 
radio flux density is highly variable and is very probably driven by 
the magnetic ac tivity whose manifestations are observed in other 
spect ral regions (jPrake et al.ll989l : lumana et al.ll993L[T995l : IWhit3 
I 2 OO 4 I) . The radio flux arises from the interaction between the stellar 
magnetic held and mi ldly relativistic particles, i.e. gyrosynchrotron 
emission (lDul]Jl985h . The same emission mechanism is at the ori¬ 
gin of radio emission from pre-main sequence (PMS) stars and 
X-ray binaries. Non-thermal radio emission also originates from 
shocks of colliding winds in massive binaries and from pulsars. 
There is growing evidence that radio flares can also occur as narrow 
band, rapid, intense and highly polarized (up to 100%) radio bursts, 
that are observed especially at low frequency (<1.5 GHz). For their 
extreme characteristics, such radio flares have generally been inter¬ 
preted as a result of coherent emission mechanisms. Coherent burst 
emission has been observ ed in different cla s ses of stellar obj ects: 


RS CVns and Elare stars JOsten et al. 

l2004ISlee et al.l2008h. Ul- 

tra Cool dwarfs JHallinan et al.l 12008 

; iRoute & Wolszczan 1201 2b 


and Chemically Peculiar (CP) stars (iTrigilio et 'Zll2008Ll2^ . All 


have, as common ingredient, a strong magnetic held, which may 
be variable, and a source of energetic particles. The number of 
stars where coherent emission has been detected is still limited 
to a few tens, because of the limited sensitivity of the available 
instruments. Thermal emission (bremsstrahlung emission) is ex¬ 
pected from winds associated with Wolf-Ray et (WR) and OB stars, 
shells surrounding planetary nebulae (PNe) and novae and je ts from 

S ’ lotic stars and class 0 pre-main-sequence (PMS) stars dWhitel 


2 THE SCORPIO PROJECT 

We have started the Stellar Continuum Originating from Radio 
Physics In Ourgalaxy (SCORPIO) project to carry out, for the first 
time, a blind deep (30 pJy) radio survey in a sky patch well suited 
for stellar work, using ATCA observations at 1.4 GHz. 

This survey has the best potential to enlarge the stellar radio 
emitting population, without suffering from selection criteria based 
on some peculiar aspects observed in other spectral bands. More¬ 
over, it will provide new insights for a better comprehension of 
the physics of particular classes of stellar systems and of plasma 
processes in a wider context. Finally, the chosen frequency (i.e. 
1.4 GHz) constitutes the best choice to answer the key question of 
how common is coherent radio emission from stellar and sub-stellar 
systems. As the space distribution of different types of stars in the 
patch selected for our survey can be assumed to be typical for the 


Galactic disk population, results from such a survey will provide us 
with a clear forecast on the potential of SKA and its pathfinders in 
the held of stellar radio astronomy. 

In addition to the scientific outcomes, results from the survey 
will be of immense value in shaping the strategy of the deep (^ 
10 pJy) surveys already planned with the SKA pathfin ders. Among 
them , the Evolutionary Map of the Universe fEMU: iNorris et aD 
l201lh . a deep, almost full sky (75 per cent) survey to be carried 
out at 1.4 GHz, is one of the major programs already approved to 
be conducted with the Australian SKA Pathfinder (ASKAP). EMU 
is expected to detect and catalogue about 70 million galaxies, and 
will also make an atlas of the Galactic Plane to an unprecedented 
sensitivity and resolution. 

The SCORPIO project will have a profound impact on the 
Galactic Plane component of the EMU survey. Specifically, it will 
guide EMU design in identifying issues arising from the complex 
continuum structure associated with the Galactic Plane and from 
the variable sources in the Galactic Plane. Moreover, it will con¬ 
tribute in evaluating the most appropriate method for source find¬ 
ing and extraction for sources embedded in the diffuse emission 
expected at low Galactic latitude. 

2.1 The selected field 

The selected field needs to satisfy the following requirements: (i) it 
should contain a sufficient number of stars, with a good spread in 
different classes of stellar objects thought to be radio emitters; (ii) it 
should contain sources already classified as radio emitters to enable 
verification of the data; (iii) it should have been already observed in 
other spectral regions to gather additional information to help with 
the classification of new, unexpected objects or to allow complete 
studies for the classified objects. 

The selected field is namely 2 x 2deg^ region centred at galac¬ 
tic coordinates / = 344.25°, b = 0.66°. Because our aim is to get 
statistical information from the survey, we choose this sky patch 
only on the a priori knowledge of the presence of a significa¬ 
tive number of stars in it. Moreover, being in the direction of the 
Scorpio constellation, it offers the possibility to probe different 
radio star populations at various distances across the Galaxy. As 
added value, in the field there is also the Sco OBI association, 
whose core is the young (^3 —5Myr) stellar cluster NGC 6231 
(jPerrv et al .111 9901) {D ~ 1.6kpc), consisting of 964 stars, where ra¬ 
dio continuum e mission from the massive sta r population has been 
already detected dSetia Gunawan et ^l2003h and X-ray properties 
of bo th OB and PMS populations have been derived (ISana et al.l 
l2006h . There are, in total, eight open clusters of stars and stel¬ 
lar associations inside the SCORPIO field, as shown in Table[T] 
[DBS2003] 176 and [DBS2003] 178 are two clusters of stars de¬ 
tected in the Two-Micron All Sky Survey (2MASS) JPutra et al.l 
l2003h . C 1653-405 and C 1658-410 are two young open clusters 
with young 0-type stars. DSH J1704.3-4204 is a cluster candidate, 
discovered in the Digiti zed Sky Survey (DSS) and 2MASS survey 
dlCronberger et al.|[2006h . that is situated in a heavily obscured part 
of the Milky Way; it could correspond to a region of lower extinc¬ 
tion rather than a physical group of star. MCM2005b 86 is a star 
cluster discovered in the Galacti c Legacy Infrared M id-Plane Sur¬ 
vey Extraordinaire (GLIMPSE) (iMercer et ^|2 oo3) . 

We have used the Simbad database to check which stellar pop¬ 
ulations are identified in the field. We found a total of 2935 objects, 
1774 of which are classified as stars (Eig.[T]). Among stars, there is 
a good spread in different classes of stellar objects thought to be 
radio emitters and a good representation of different spectral types. 
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Table 1. Stellar Clusters and Associations in the SCORPIO field 


Name 

RA 

(J2000) 

Dec 

(J2000) 

1 

(deg) 

b 

(deg) 

NGC 6231 

16 54 08.5 

-41 49 36 

343.4601 

-rOl.1866 

Sco OB 1 

16 53 28.8 

-41 57 00 

343.2900 

+01.2000 

[DBS2003] 176 

16 59 23.0 

-42 34 23 

343.4830 

-00.0380 

[DBS2003] 178 

17 02 10.1 

-41 46 48 

344.4260 

+00.0440 

C 1653-405 

16 57 00.0 

-40 40 01 

344.7010 

+01.4960 

C 1658-410 

17 02 12.0 

-41 04 12 

344.9900 

+00.4800 

DSH J1704.3-4204 

17 04 20.1 

-42 04 24 

344.4394 

-00.4546 

MCM2005b 86 

17 04 40.1 

-41 42 25 

344.7680 

-00.2810 



Spectral Type 



Clouds 

Stars 

1 1 1 1 1 111 1 1 1 1 1 1 11 i I 1 1111 l,.J 1 1 I 

1 10 100 1000 
Number of Objects 

Figure 1. Histogram showing the number of each different possible radio 
emitting class, listed in SIMBAD, falling in the selected field. The labels 
follow the SIMBAD nomenclature. There are many sources with only gen¬ 
eral X-ray, UV and IR classification. Bubbles, SNR, y-ray sources. Pulsars, 
PN, H II regions and High and Low mass X-ray binaries are also repre¬ 
sented. For 14 sources we found only generic radio source classification 
(radio, mm and sub-mm sources). 


Part of the prop osed patch has already been surveyed in the 
mid-IR by SPITZER (iBeniamin et al.ll2003l:JCarev et al.lf2009h and 
by Herschel HI-GAL survey JlVIolinan et al.l20icl) and will be cov¬ 
ered in the near future by CORNISI^^ south, that will nicely com¬ 
plement the radio spectral information for the brightest detections. 


^ CORNISH south is a high resolution survey of the Galactic Plane (PI: M. 
Hoare) being carried out with ATCA at 6 and 9 GHz with a target rms noise 
of 0.18mJybeam“^ 


Table 2. Log of Observations 


Date 

Configuration 

time on source 



(h) 

2011 April 21-24 

6A 

38 

2012 June 3-11 

6B 

6 


Figure 2. Distribution of the subsample with Spectral Classification as re¬ 
ported in Simbad as function of spectral types. The larger number of B-type 
stars is a consequence of the presence of a part of the ScoOB 1 association 
and of its nuclear cluster NGC 6231 in the field. 

3 OBSERVATIONS AND DATA REDUCTION 
3.1 The pilot experiment 

A pilot experiment for the SCORPIO project took place in 2011, 
from April 21 to 24 and in 2012, from June 3 to 11 (see Tab[2]). 
In this experiment, only about one quarter of the selected field was 
observed, namely a 0.5 x 2deg^ region centred at I = 343.5°, b = 
0.66°. The field was observed with ATCA, in mosaic mode, using a 
8.8 arcmin spacing hexagonal grid, requiring a total of 38 pointings. 

We cycled the array around these pointing centres, spending 
about 1 min at each pointing centre. The secondary calibrator 1714- 
397 was visited every 20 min to calibrate the complex antenna 
gains. PKS B1934-638, which is the standard primary calibrator 
for ATCA observations, was observed at the beginning of each ob¬ 
serving run. Each pointing was observed for 1.16h in total. The 
total observing time for all pointings was 58 h, including a 20 per 
cent overhead for fiux, bandpass and phase calibration. 

The observations were made with the new broadband back¬ 
end syste m for the ATCA, th e Compact Array Broadband Backend 
(CABB) dWilson et al.ll201lh . with an effective observing band of 
2 GHz divided into 2048 1-MHz channels. All four Stokes parame¬ 
ters were measured. 

REI strongly corrupted our data, mostly between 1.2 and 
1.4 GHz, requiring a large amount of flagging at those frequen¬ 
cies. Elagging and calibration were performed in MIRIAD, using 
the MIRELAG task. We ended up flagging about 30-40 per cent 
of the data. The usable range of frequencies goes from 1.350 to 
3.100GHz, corresponding to a frequency variation of about 80 per 
cent across the band. 

Calibration was performed by using the standard MIRIAD 
tasks. The well studied Gigahertz Peaked Spectrum radio galaxy 
PKS 1934-638 was used as b andpass and flu x calibrator. A flux 
density of 12.31 Jy at 2.1 GHz (iRevnoldsll 19941) was used to derive 
the flux of the phase calibrator J1714-397 (2.031 Jy with an accu¬ 
racy of 2%). 


3.2 Map making 

Imaging was performed in MIRIAD and the deconvolution was 
done with the task MECLEAN using the Hogbom algorithm 
(lHogbomlll974l) . All data for each pointing were combined into 
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Figure 3. Final map of the pilot field of Scorpio showing the ATCA field of view for each pointing at 1.470GHz (left panel) and 2.930GHz (right panel). 
The region chosen for source finding, including all the frequency sub-bands, is delimited by the red line. Grid coordinates are Galactic. Directions are RA 
(horizontal)-Dec (vertical). 


Table 3. Frequency sub-bands and primary beam radii at 5% of the maxi¬ 
mum amplitude 


sub-band 

central v 
(GHz) 

Av 

(GHz) 

R(5%) 

(arcmin) 

1 

1.469 

0.169 

30.8 

2 

1.649 

0.190 

27.5 

3 

1.850 

0.213 

24.4 

4 

2.075 

0.239 

23.6 

5 

2.329 

0.268 

20.1 

6 

2.613 

0.300 

18.7 

7 

2.932 

0.337 

16.7 


one file. Each pointing was then imaged separately prior to mo¬ 
saicking. Five iterations of self-calibration were applied, the first 
three with phase calibration only and the last two with both phase 
and amplitude calibration, progressively increasing the number 
of CLEAN components used to model the sky emission. Multi¬ 
frequency CLEAN was used to model the spectral variation of the 
source emission using a two-term Taylor-Polynomial. The uv data 
were weighted by adopting a robust weighting, with parameter 0.5. 
This permitted us to reach a better signal-to-noise ratio at the ex¬ 
pense of a poorer angular resolution when compared to a uniform 
weighting scheme. 

Next, the self-calibrated uv data were divided into seven sub¬ 
bands, each of the same fractional bandwidth Av/v = 10%, im¬ 
aged separately and adjusted to the reference position a = 16:56: 
28.912, 5 = —42°03'31.28". Splitting up the data into seven sub¬ 
bands serves to improve the accuracy of the primary beam correc¬ 
tion. In TableO the central frequency, bandwidth and primary beam 
size for each sub-band are listed; R(5%) is the radius of the primary 
beam out to a level of 5% of the maximum, following the Gaussian 


prima ry beam model for ATCA telescopes (Iwieringa & Keste\^ 

The uv data were tapered in each sub-band differently to ob¬ 
tain nearly identical synthesized beams. The cell size was set to 
1.5arcsec, suitable for a good sampling of the synthesized beam. 
Pointings were CLEANed to a depth of 5 a and restored using a 
Gaussian beam of size 14.0 by 6.5 arcsec. 

A mosaic for each sub-band was formed using the LINMOS 
task. Pointings were corrected for the primary beam and mosaicked 
together, weighting them by their respective rms noise values. The 
bw option was used to average the primary beam response over the 
sub-band (rather than use the primary beam response at the central 
frequency of the sub-band). Finally, the seven sub-band mosaics 
were combined using the IMCOMB task, assigning equal weights 
to the sub-bands. 

In Fig. 0 the final map of the pilot experiment is shown. It 
is overlayed with the field of view (defined as the region out to 
5% of the primary beam maximum) of the ATCA antennas and 
centred on the pointing positions used in the mosaiced observations 
(see Tab. O. In the left panel the field of view refers to the lowest 
frequency sub-band (v ^ 1.470 GHz) and in the right panel to the 
highest frequency sub-band (v 2.930GHz). 


3.3 Mapping extended sources 

The standard CLEAN algorithm is an iterative procedure which as¬ 
sumes that the sky is a collection of point sources on an empty 
background. The algorithm looks for all these sources and creates 
a ‘CLEAN components’ list, which is used to obtain a reliable rep¬ 
resentation of the sky by means of a convolution with an ideal PSF. 
This procedure models an extended source as a set of point sources. 
However the results may not be completely satisfactory. 

Several algorithms have been proposed to properly model ex¬ 
tended sources also. These methods assume that the sky is not com- 
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Figure 4. The trend of rms (from the noise map) as a function of the Galac¬ 
tic latitude for both Stokes / and V (note the use of a logarithmic scale for 
the j axis). For each pixel of the map the rms has been averaged over 0.3 deg 
in longitude. The increase in the rms close to the Galactic Plane in Stokes I 
indicates the presence of diffuse emission. 

posed of point sources only, but of sources of many different sizes 
and scales. 

In pa rticular, for our p urpose, we used the multi-scale CLEAN 
algorithm dCornwe lll2008h . as implemented in CASA. With respect 
to the usual CLEAN procedure, we defined three different scales to 
model source extension: 0 to safely search for point sources (equiv¬ 
alent to a normal CLEAN), 5 pixels (about the beam dimension) 
and 15 pixels. A few examples of extended sources, imaged us¬ 
ing the multi-scale CLEAN, are presented in Eigs. [5]and[6| Multi¬ 
scale CLEAN was found to be much more effective than standard 
CLEAN at removing the sidelobes around these sources. 

4 IMAGE ANALYSIS AND SOURCE CATALOGUE 

4.1 Source finding 

Sources were detected through an analysis of the mosaic combin¬ 
ing all the frequency sub-bands. A significant non-zero background 
is expected if there is large-scale, diffuse emission, which is the 
case in the SCORPIO field given its location on the Galactic Plane. 
Since our aim was to detect and characterize compact sources in the 
field, the background was measured and subtracted from the image 
before source extraction. The background was estimated for each 
pixel using the following procedure: because computing the back¬ 
ground is fairly computationally expensive, it was only calculated 
for a limited number of pixels occuring at regular intervals in the 
V and y directions. Eor each of these pixels, the background was 
taken as the median pixel value inside a square centred on the pixel 
with a width of 10 times the synthesized beam, but a standard itera¬ 
tive technique was used where points were clipped if they deviated 
from the median by more than three times the rms of the devia¬ 
tion from the median in the previous iteration. The process was 
terminated when the rms of the deviation from the median from 
successive iterations differed by less than 10%. The background in 
the remaining pixels was estimated using a bicubic interpolation. 

Having subtracted the background from the image, the next 
step involved determination of the noise level. Eor each pixel, an 
initial estimate of the noise was taken as the rms inside a box of 


size 10 times the synthesized beam, centred on the pixel. To pre¬ 
vent the noise estimate from being affected by real source emis¬ 
sion, pixels outside the range d=3(7 were removed and the rms was 
re-calculated. This process was repeated a number of times until 
the noise was found to decrease by no more than 10 per cent. 

Sources were detecte d and characterized using a similar 
method to that employed bv lEranzen et al.l ([2oTl|) for the IOC sur¬ 
vey. The noise maps were used to identify sources on the basis of 
their signal-to-noise ratio. Local maxima above 5a were identified 
as sources. A peak position and fiux density value were measured 
by interpolating between the pixels. A centroid position, integrated 
flux density and source area were also calculated by integrating 
contiguous pixels down to 2.5(7, and sources were identified as 
overlapping if the integration area contained more than one source. 

4.2 Sensitivity 

There are a number of factors which affect the sensitivity of our 
SCORPIO mosaic: the effective bandwidth, system temperature, 
integration time and other instrumental characteristics, as well as 
the ability to deconvolve extended sources, which in turn depends 
on the uv plane sampling, in particular on the minimum baseline of 
the interferometer. The Galactic Plane is full of extended sources 
and diffuse synchrotron emission. In order to quantify the theoret¬ 
ical sensitivity limit, we made a V polarization map of the whole 
field with the aim to compare the Stokes I and V rms maps. 

In general, radio source emission presents a circular polariza¬ 
tion close to 0%. We can assume that no source is present in the V 
map which therefore does not contain artifacts due to deconvolu¬ 
tion errors and non-perfect uv sampling. Stokes I and V are given 
by linear combinations of the cross correlations of the responses of 
the feeds, which are linearly polarized for ATCA {X vertical and Y 
horizontal polarizations): 

2I = XX^YY (1) 

and 

2iV=XY-YX. (2) 

Since all the cross products have the same statistics, I and V maps 
have the same statistics too, and therefore the noise level should be 
comparable. Pig. |4] shows the rms for the I and V maps. At Z? ^ 0° 
the I map has a very high rms compared with the V one. This is 
imputed to the presence of extended sources and diffuse emission in 
the Galactic Plane not correctly sampled in the wv-plane due to the 
lack of short baselines. This leads to fiuctuations which adds noise 
into the map. At the edge of the field (b < —0.25° and b > 1.8°) 
the noise increases due to the primary beam attenuation and the 
effective frequency of the image decreases due to the frequency 
dependence of the primary beam (see Pig. [3]). The rms of the V 
map is in general lower, reaching about 20pJybeam“^, while for 
the I map it does not go below 30 pJybeam“^. In any case the high 
brightness of extended synchrotron emission at \b\ < 0.4° {Tb ~ 
10 K at 16 cm) increases the rms of any polarization map. 

4.3 The compact components catalogue and Cross 
Identifications 

Radio components were extracted from the image as described in 
subsection l4.1l Sources were extracted only from a rectangular sub- 
region of SCORPIO with the best sensitivity. This is the region that 
includes the high frequency pointings and is defined by the rectan¬ 
gle between I = 343° and I = 344° and b = —0.5° and b = 1.833°, 
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indicated in red in Fig.|3l Components corresponding to or affected 
by artifacts were removed. The resulting catalogue contains 614 
radio components; Table |4] lists the first 25 components in the cat¬ 
alogue while the complete catalogue is available in the online ver¬ 
sion. Querying SIMBAD and NED (search radius of 6 arcsec) we 
found 34 matches. Some of these are classified as stars, PN or 
Pulsars, others have only a generic classification (i.e. radio or IR 
source). 

Among those 34 matches, we individuate 10 sources that 
are associated with objects that could be stars, 2 Planetary Neb¬ 
ulae (SC0RPI01_614 and SC0RPI01_325) and 2 Pulsars (SCOR- 
PI01_399 and SCORPIOl_091). In the following list, more infor¬ 
mation on the possible stars are provided. The classification of 
these objects, as reported in the SIMBAD database, is indicated 
between square brackets. 

[*] SC0RPI01_534: generically classified in SIMBAD as star, 
no references are reported in the SIMBAD database. The source 
SCORPIO 1_534 is embedded within extended radio emission, (see 
Section [5^ . 

[IR*] SCORPIO 1_012: associated to Sco, classified by 
SIMBAD as Variable Star of irregular type. 

[WR*] SCORPIO 1_118: WR star, cha racterized by strong 
stellar wind already known as radio source (iBieging et al.lll982l : 
lAbbott et ahlll9861 : feloggll 198^ . The w ind radio emission has been 
detected at 1.3 mm (11.6 ±2.3 mJy) (iLeitheree et al.|[l9^ . The 
spectral index, estimated between our SCORPIO measurement and 
the millimeter measurement, is consistent with radio emission from 
a stellar wind. 

[*] SC0RPI01_311 (spectral type: WC7-I-06V): well known 
WR st ar (HR 6265/HD IS^VOAVR 79) already known a s radio 
source (IBieging et aPl 19821:lAbbott et al.lll98fil:lHogglll989h . 

[iC*] SCORPI01_313 (Cl* NGC 6231 SBL 489): classified 
as Star in Cluster, and was list ed in the photoin etric study of the 
young open cluster NGC 6231 (ISung et al.lll998h . 

[Y*?] SC0RPI01_468 (G343.7018-h00.086): classified by 
SIMBAD as Young Stellar Object Candidate. 

[*] SC0RPI01_243 (IRAS 16495-4140): classified as star in 
the SIMBAD database, no references reported. 

[Be*] SC0RPI01_219 (V921Sco): well known Herbig Ae/Be 
star (128 references listed by SIMBAD). This star is surrounded 
by an extended envelop detected at the sub-mm wavelengths 
(lManningslll995h and is also visible in the infra red. 

[*iC] SCORPIOl_406 (Cl* NGC 6231 SBL 759 ): classified 
as Star in Cluster, li sted in the photom etric study of the young open 
cluster NGC 6231 (ISung et al.lll998h . The compact source SCOR¬ 
PIO 1_406 is embedded within more extended radio emission. 

[*iC] SCORPIOl_350 (Cl* NGC 6231 BVF 77: (spectral 
type: B8.5V): classified as Star in Cluster and listed in the pho¬ 
tomet ric study of the young open cluster NGC 6231 dSung et al.l 
Il998h . 


The first step in validating and classifying the detected radio 
sources, with no matches in SIMBAD or NED, is to determine how 
many of these have counterparts at other wavelengths. Hence, we 
cross-matched our radio sources with major catalogues, from the 
optical to the radio frequency range, covering the SCORPIO re¬ 
gion, name ly: the Naval Observat ory Merged Astrometric Dataset 
(^NOMAD: IZacharia s^t^ 2004 ). the 2MASS all sky catalogue 
of point sources dCutri et al .1120031) . the Galactic Legacy Infrared 
Mid-Plane Survey Extraord inaire Source Catalogue (GIMPSE; 
ISpitzer Science Centei]l2009h. the Wide-field Infrared Survey Ex¬ 
plorer (WISE; Cutri et all l2014h . the Midcourse Space Experi- 


Table 4. Number of matched sources 


Catalog Name 

Number of Matches 

NOMAD 

320 

2MASS 

301 

GLIMPSE 

229 

WISE 

116 

MSX6C 

47 

AKARI 

34 

Hi-Gal 

148 

IRAS 

117 

ATLASGAL 

14 

MGPS-2 

43 

WBH2005 

18 

RMS 

6 


ment infrared point source catalogu e (MSX6C; lEgan et al.ll2003h . 
the AKARI mid-IR (AKARI/IRC; llshihara et al.l l201Qh all sky 
survey, th e Herschel infrared Galactic plane Survey catalogue 
(Hi-GAL; iMolinari et al.ll201Qll . the IRAS point source catalogue 
(IHelou & Walked 19881) . the APEX Telescope L arge Area Sur- 

S if the GALaxy at 345 GHz (ATLASGAL; IContreras et al.l 
, the 2nd Epo ch Molonglo Galactic Plane Surveys (MGPS-2; 
Murphy et al .l2007h . the cat alogues of compac t radio sources in the 
Galactic Plane (WBH2005;|m te et al.l2005l) and th e radio survey 
of the Red MSX Source (RMS; Urquhart et al.ll2007h . 

In performing such a task the most critical information is posi¬ 
tion accuracy. The position errors in right ascension and declination 
for the radio s ource s were derived following equations 3, 4 and 5 in 
IPurcell et~^ d2013l) . that take into account the effective signal-to- 
noise ratio for each point source and the interferometer beam size, 
assuming for the sky background a Gaussian correlated noise. 

An error of 0.1 arcsec, co rresponding to the position uncer¬ 
tainty of the phase calibrator dNorris et al.ll2006h . was added in 
quadrature to this value, to obtain the final error on the radio posi¬ 
tion (cTradio)- the cross-correlation process, we set a search radius 

of r = Sdposi, where (jposi = and CTcat is the position 

error in the reference catalogues (and all other instances in paper). 

The number of sources in our catalogue that match the entries 
in the reference catalogues is summarized in Table |4] Among the 
614 SCORPIO sources, 487 have a counterpart in at least one of the 
inspected catalogues. There are no SCORPIO sources concurrently 
present in all inspected catalogues, however: 78 per cent of SCOR¬ 
PIO sources found in one of the IR catalogues (GLIMPSE, WISE, 
MSX) are also in NOMAD or 2MASS; 53 per cent of those found 
in one of the far-IR catalogues (AKARI, IRAS, Hi-GAL, ATLAS¬ 
GAL) are also in NOMAD or 2MASS; only 34 per cent of those 
found in one of the radio catalogues (RMS, WBH2005, MGPS-2) 
are also in NOMAD or 2MASS. 

The Compact Components Catalogue (Table [5]) is organized as 
follows: 

Column (1) Component number. This is the internal definition of 
the component used within this paper. If two or more sources have 
the same number (with different letters) they are components of a 
group, i.e. associated with the same radio source, as returned by 
the adopted source finding method (see Section 4.1) and confirmed 
by a visual inspection. 

Column (2) Source Name; 

Column (3) Galactic longitude of the component; 

Column (4) Galactic latitude of the component; 
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Figure 5. Composite picture of the field centred on [DBS2003] 176. The 
sub-panel shows only the SCORPIO map, while in the background panel the 
mid-IR/FIR maps from Spitzer (IRAC, 8 p.m, green) and Herschel (PACS, 
70 p.m, red) are superimposed on the SCORPIO map, in blue. The white 
arrow indicates the position of the the compact component SCORPIO 1_300. 
The white circle points out that there is no radio emission associated with 
S16 (see text). 

Columns (5) and (6) Right ascension and Declination (J2000.0) of 
the peak of the emission; 

Columns (7) and (8) integrated flux density at 20 cm 

(S') and its estimated uncertainty ( AS), defi ned _ as 

y^rms^ + (S x AW/W)^ + in mJy (iFomalontI Il989h : 

where rms is the local rms, W and AW are the beam width and 
its associated error (AW = rms xW/P, P is the peak intensity) 
and a is the relative error of the flux density of the calibrator. We 
assumed a = 0.03; 

Column (9) Area of the component in beam units; 1 means source 
not resolved; 

Column (10) match with source in other catalogue. N: NOMAD, 
M: 2mass, G: GLIMPSE, W: WISE, X: MSX, A: AKARI, H: 
HIGAL, I: IRAS, a: Atlasgal, m: MGPS-2, w: WBH05, R: RMS. 


5 EXTENDED SOURCES: A FEW EXAMPLES 

Our source extraction method is biased against extended sources. 
However, a visual inspection of the flnal map points out the pres¬ 
ence of 17 of this kind of sources falling in the SCORPIO held. In 
particular, there are many round ish structures recalling what sev¬ 
eral authors named as ’’bubbles” dChurchwell et al.ll2006h . Bubbles 
are pervasive throughout the entire Galactic Plane and have been 
mainly discovered by means of t he new extended surve y conducted 
with Spitzer, namely GL IMPSE (iBeniamin et al.ll2003h and MIPS- 
GAL dCarev et al.l2009l> as they are quite often associated with ex¬ 
tended dusty structures. 

Infrared bubbles are usually associated with hot young stars 
in massive star formation regions and in some cases they are co¬ 
incident with known H II regions (IChurchwell et al.ll2006l) . How¬ 
ever, there is growing observational evidence that some bubbles 
can be related to massive evolved stars, such as LB V, WRs and 
SNRs (IWachter et aP 120101 : iGvaramadze et al.l 120101) or with PNe 
(lingallinera et al.ll20l4) . Despite the apparent similar morpholo¬ 
gies, there are many differences between the emission character¬ 
istics which stand out clearly when different maps, obtained at dif¬ 


ferent IR and radio bands, are compared. As an example, bubbles 
associated with massive star formation regions have an extended 
8 qm emitting region embracing the more compact 24 qm one, usu¬ 
ally co-spatial with the ionized region, traced by the radio. This is 
interpret ed as a signature of a PD R, traced by the PAHs emission 
at 8|am (IChurchwell et al.ir2006h . In contrast, the emission char¬ 
acteristic of a Bubble associated with an evolved star has a more 
complex structure. The dusty envelope i s usually brighter and more 
extended at 24 qm (IWachter et al.ll201Ql) but the ionized region can 
be well contained or c o-spatial or even more extended with respect 
to the 24 \im emission jChu et al.ll2009h . This is interpreted as a re¬ 
sult of the very complex dust distribution in the circumstellar enve¬ 
lope as consequence of different mass-loss episodes that can occur 
during the evolution of the central star. 

In the following we will provide a few examples of extended 
sources detected in the SCORPIO pilot. Some of these are associ¬ 
ated with already known Galactic sources, while others are unclas- 
sifled. Our work provides the first high resolution radio continuum 
maps for most of these sources, revealing the details of the ionized 
gas. 

In presenting our results we also consider the dusty environ¬ 
ment associated with the radio sources and display the radio maps 
togeth er with maps of the r egions provided by GLIMPSE and HI¬ 
GAL (iMolinari et al .11201^ . that trace the warm and cooler dust 
distribution in the surrounding regions. 

5.1 [DBS2003] 176 

The stellar cluster [DBS2003] 176 is associated with a mul- 
tiple bubble (S16 and S17) detected in the GLIMPSE survey 
dChurchwell et all l2006l) . The H II region associated with the 
bubble S17 has been detec ted in the Parkes-MIT-NRAO (PMN) 
survey dWriehtet all 11994 at 4.5 GHz, which has a beam size 
of 4.1arcmin. Radio emission at 20cm from an extended non- 
Gauss ian source has also been reported bv IZoonematkermani et al.l 
dl990h . 

Our SCORPIO map reveals, for the first time, the extended 
radio emission at 20 cm, with unprecedented detail, thanks to an 
angular resolution comparable to mid-IR/FIR images. In Fig [5] the 
radio emission (SCORPIO) is pictured in blue, the Herschel (Hi- 
GAL) map at 70 qm in red, and the Spitzer (Glimpse) map at 8 qm 
in green. The extended radio emission is well contained inside the 
dust emission, traced by the IR. In the false color image the asso¬ 
ciated stellar cluster is recognizable inside bubble S17. This stellar 
cluster is spatially resolved in the mid-IR (8 qm) in many stellar 
components. 

The two bubbles, S16 and S17, are clearly visible in the 
IR maps, while there is no extended radio emission associated 
with bubble S16 (sub-panel). Two point-like radio sources (SCOR¬ 
PIO 1_320 and SCORPIO 1_300) have been retrieved by our extrac¬ 
tion procedure. The position of SCORPIO 1_300 in the compos¬ 
ite RGB map is indicated by the arrow, see FigO The brighter 
radio point source (SCORPIO 1_320), clearly visible within the 
two IR bubbles S16 and SI7, is an already known radio source 
JZoonematkermani et al.l ll99'Qh and is present in The Two Mi- 
cron All Sky Survey (2MASS) catalogue of extended sources 
dSkrutskie et al.l 1200^ . A flux density measurement at 5 GHz 
(176.9 mJy) has also been reported for SCORPIO 1_320 and the 
source has been tentatively classified as a Massi ve Young Stellar 
Object (MYSO) candidate dUrauhart et al.ll200^ . Source SCOR¬ 
PIO 1_300 is located within bu bble S17 and has already been de- 
tected at radio wavelengths by IZoonematkermani et al.l (ll99C)h . It 
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Table 5. First 25 entries of the Point Source Catalogue of the Pilot experiment. The complete catalogue is available in the online version. 


ID 

1 

(deg) 

b 

(deg) 

RA 

(J2000) 

Dec 

(J2000) 

(mJy) 

Ay 

(mJy) 

Area 

(beam) 

Matching 

SCORPIOLOOl 

343.0025 

1.7604 

16:50:12.21 

-41:48:56.2 

33.33 

1.03 

1.0 

- 

SCORPIOl_002 

343.0051 

0.2234 

16:56:39.13 

-42:47:08.2 

0.49 

0.09 

1.0 

- 

SCORPIOl_003 

343.0134 

0.6086 

16:55:02.85 

-42:32:15.5 

7.74 

0.29 

1.4 

NMGm 

SCORPIO 1_004 

343.0138 

1.7208 

16:50:24.33 

-41:49:56.4 

1.60 

0.17 

1.0 

N 

SCORPIOl_005a 

343.0139 

1.1508 

16:52:46.37 

-42:11:42.1 

12.80 

0.50 

2.0 

N 

SCORPIOl_006 

343.0152 

0.1166 

16:57:08.49 

-42:50:39.6 

1.27 

0.16 

1.7 

NMGW 

SCORPIO 1_007 

343.0157 

-0.1830 

16:58:25.48 

-43:01:49.4 

1.51 

0.12 

1.2 

NMG 

SCORPIOl_005b 

343.0186 

1.1577 

16:52:45.61 

-42:11:13.0 

106.56 

3.20 

1.4 

- 

SCORPIO 1_009 

343.0201 

1.5248 

16:51:14.28 

-41:57:09.1 

1.34 

0.11 

1.0 

- 

SCORPIOLOlO 

343.0203 

0.3109 

16:56:19.93 

-42:43:08.3 

20.97 

0.64 

1.2 

Gmw 

SCORPIOLOll 

343.0216 

0.7407 

16:54:31.13 

-42:26:53.1 

1.12 

0.10 

1.3 

G 

SCORPIOl_012 

343.0277 

0.8702 

16:53:59.71 

-42:21:42.3 

1.17 

0.09 

1.2 

NMXAH 

SCORPIOLOIS 

343.0317 

0.8579 

16:54:03.63 

-42:21:59.3 

0.44 

0.08 

1.0 

NMH 

SCORPIOl_014 

343.0322 

0.9655 

16:53:36.65 

-42:17:53.0 

1.21 

0.12 

1.3 

NH 

SCORPIOl_015 

343.0356 

0.2182 

16:56:46.68 

-42:45:53.7 

8.06 

0.25 

1.0 

N 

SCORPIOl_016 

343.0394 

1.5983 

16:50:59.99 

-41:53:26.8 

2.49 

0.13 

1.0 

I 

SCORPIOl_017 

343.0415 

1.2316 

16:52:31.80 

-42:07:20.7 

1.46 

0.13 

1.0 

NM 

SCORPIOLOIS 

343.0428 

0.0266 

16:57:37.15 

-42:52:43.9 

5.29 

0.22 

1.0 

- 

SCORPIOl_019 

343.0457 

0.8848 

16:53:59.73 

-42:20:18.9 

71.55 

2.15 

1.0 

Hm 

SCORPIO 1_020 

343.0461 

1.2752 

16:52:21.85 

-42:05:28.5 

4.71 

0.25 

2.2 

I 

SCORPIOl_021 

343.0468 

0.3010 

16:56:27.87 

-42:42:16.1 

0.47 

0.06 

1.0 

NMGW 

SCORPIOl_022a 

343.0489 

0.9666 

16:53:39.77 

-42:17:04.1 

109.20 

3.28 

1.6 

Gm 

SCORPIOl_023 

343.0501 

1.6133 

16:50:58.45 

-41:52:22.7 

1.75 

0.17 

2.1 

NMWAI 

SCORPIOl_024 

343.0502 

0.4249 

16:55:57.01 

-42:37:27.4 

1.04 

0.09 

1.0 

NMGH 

SCORPIOl_022b 

343.0514 

0.9634 

16:53:41.09 

-42:17:04.3 

23.35 

0.72 

1.6 

m 



Figure 6. Composite picture of the held around IRAS 16566-4204. In 
the held two radio sources are visible. The point-like source is associ¬ 
ated with SCORPIO 1_600, while the more extended source is associated 
with SCORPIO l_592a and SCORPIO I _592b. The mid-IR/FIR maps from 
Spitzer (IRAC, 8 p.m, green) and Herschel (PACS, 70p.m, red) are superim¬ 
posed to the SCORPIO map (blue) with colour coding as indicated. 


is also associated with a GLIMPSE source. From Fig|5] SCOR¬ 
PIO 1_300 seems to be loeated in the top region of a dust pillar, well 
visible at the IR wavelengths, indicating possible triggered star for¬ 
mation by the expanding H II region. Further investigation of that 
source is necessary to identify it as a YSO. 


5.2 IRAS 16566-4204 

In Fig.[6l the held around the position of IRAS 16566-4204 is 
shown. Two radio sources are clearly visible, one very compact 
(SCORPIO 1_600) and the other more extended. Their positions are 
eoineident with those of IRAS 16566-4204. IRAS 16566-4204 has 
been previousl y classided as a UCHII on th e basis of IRAS colour- 
colour criteria (IWood & Churchw^ll989h an d of the associatio n 
with both radio continuum and maser emission dWalsh et al.l 19971) . 
However, higher resolution ATCA observations revealed a compact 
radio source and a masing methanol (6.668 GHz) sour ce with a po¬ 
sition offset with respeet to the IRAS pointing center dWalsh et al.l 
Il998h . Such a picture together with the SCORPIO map is consistent 
with the UCHII being associated with SCORPIO 1_600 and IRAS 
16566-4204 associated with the extended radio source whose elas- 
sihcation is not clear. However, the radio emission is well contained 
inside the more extended 8 p.m emission (see Fig[6]), providing us 
with a morphological hint to classify IRAS 16566-4204 as a mas¬ 
sive star formation region. At the position of IRAS 16566-4204 a 
101 mJy souree at 4.85 GHz is reported in the Parkes-MIT-NRAO 
(PMN) survey. 


5.3 IRAS 16573-4214 


The souree has been elassihed as a bona hde high-mass protostel- 
lar object on the basis of its IR colo urs. The source was observed 
at 6.7 GHz bv iMacLeod et al.l (Il998h to search for methanol maser 
but this was not detected. The sou rce was instead dete eted in the 
molecular lines of CS and C^^O bv iFontani et al. (l2005h. iinplying 
its association with molecular clouds. ISanchez-Monge et al.] J2OI3I) 
reeently observed the field with ATCA to seareh for H2O maser 
emission. While radio continuum emission at 18 and 22.8 GHz was 
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Figure 7. Composite picture of the field centred at the position of 
IRAS 16573-4214. The mid-IR/FIR maps from Spitzer (IRAC, 8 p.m, green) 
and Herschel (PACS, 70p.m, red) are superimposed on the SCORPIO map 
(blue) with colour coding as indicated. 



Figure 8. Composite picture of the field centred at the position of IRAS 
16561-4207. The mid-IR/FIR maps from Spitzer (IRAC, 8p.m) and Her¬ 
schel (PACS, 70p.m) are superimposed on the SCORPIO map with colour 
coding as indicated. 


detected at the IRAS position, an H 2 O maser was detected 20" 
to the edge of the continuum emission. Continuum emission at 
1.2 mm co incident with the H 9 O maser was detected with the SEST 
telescope (iBeltran et al.ll2006h . 

The field around the IRAS source IRAS 16573-4214 is shown 
in Fig. |7] as a superimposition of the radio image on mid-IR/FIR 
maps from Spitzer (IRAC, 8 p.m) and Herschel (PACS, 70 p.m). The 
SCORPIO map, obtained at 2.1 GHz with an improved spatial res¬ 
olution with respect to those bv ISanchez-Monge et al.l (l2013l) . re¬ 
veals much more detail of the radio source morphology and in par¬ 
ticular a disk-like shape, well embedded in a dusty environment as 
traced by the 70 jam emission. The H 2 O maser emission is not asso¬ 
ciated with IRAS 16573-4214 and we do no t detect any radio con¬ 
tinuum emission at the position reported bv ISanchez-Monge et akl 
il2013h . at our local sensitivity threshold (5 cr) of 500 pJy. 

5.4 IRAS 16561-4207 

In Fig| 8 ] the SCORPIO region associated with IRAS 16561-4207 
is shown superimposed on mid-IR/FIR maps from Spitzer (IRAC, 
8 jam) and Herschel (PACS, 70p.m), with colour coding as indi¬ 
cated. The brightest source on the left is SCORPIO 1_534, classified 
as a generic star in the SIMBAD database with no references asso¬ 
ciated. The source is imbedded in extended radio and IR emission. 
Both 8 p.m and 70 jam share the same morphology as well as the 
same extension, while the radi o emission appears to be well con¬ 
tained inside the dusty cocoon. IWalsh et ^ J2014h reported the de¬ 
tection of aH 20 maser spot at a = 16:59:34.3, S = — 42°ll'59.8". 
It is not clear if this maser spot is associated or not to the source, 
being localized in the external border of the IR emission. 

5.5 IRAS 16520-4146 

In FigH the SCORPIO field around IRAS 16520-4146 is shown 
superimposed on mid-IR/FIR maps from Spitzer (IRAC, 8 jam) and 
Herschel (PACS, 70p.m), with colour coding as indicated. To bet¬ 
ter appreciate the radio morphology of the sources detected in the 
SCORPIO field, only the radio data are shown in the insert in the 
bottom right corner of the same figure. 


In the SCORPIO map two extended components are evident: 
the one in the south-west is associated with IRAS 16520-4146, 
while the source in the north-east is not classified in the SIMBAD 
database but is reported in the Mizuno catalog of mid-IR Galactic 
bubbles. Moreover, while IRAS 16520-446 has a mid-IR/FIR coun¬ 
terpart, this is not the case for the other source in the field, that how¬ 
ever shows a clear roundish morphology at 24 p.m (iMizuno et al.l 
ImI. Interestingly, a hard X-ray transient was detected by INTE¬ 
GRAL (IGR J16 558-4150) at a = 253.95°, S = -41.83°, with 3- 
arcmin accuracy (ISoldi et al.ll2007h . No radio counterpart is visible 
at the coordinates of the INTEGRAL source. However, due to the 
low accuracy of the INTEGRAL position we cannot exclude that 
the unclassified radio source, detected to the north-east of IRAS 
16520-446, is related to the X-ray source. Further investigation is 
necessary. 


6 SOME PRELIMINARY CONCLUSIONS 

The aim of this paper is to provide an overview of the SCORPIO 
project and to present the results of the pilot observations. Some 
conclusions can however be drawn on the basis of these prelimi¬ 
nary results. In spite of the large number of objects located in the se¬ 
lected field and classified as stars, according to SIMBAD, we found 
very few matches that definitively associate SCORPIO radio com¬ 
ponents with stellar objects. Among the matches with the SIMBAD 
database, 10 sources are classified as stars. All of them appear re¬ 
lated to thermal emission and associated with an early-type star (see 
Section 1431 ). However, this could be an effect of the higher number 
of stars with spectral classification ‘B-type’ in the selected field as 
indicated by the distribution of spectral types shown in Fig.|2] 

The mere presence in the field of stars belonging to a class 
thought to be radio emitting is a necessary but not sufficient condi¬ 
tion to detect them, as a series of requirements needs to be satisfied: 
a sufficient mass-loss rate and ionizing UV field, for thermal emit¬ 
ters, and a sufficient magnetic field strength and energetic particle 
density, for non-thermal ones. Distances also play their role and 
variability is expected in most of the non-thermal radio emitters. 

All the above considerations translate into a real difficulty to 
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quantify the number of stellar objects we expect to detect in a spe¬ 
cific patch of the sky. 

We have detect ion rates going from 20 per cen t for OB stars 
teieg ing et al.lll989l) and 25 per cent for CP stars (iTOgdio^et^ 

g to 30-40 per cent for the active binary systems dUmana et al. 

. However, all information collected so far comes from small 
samples biased towards some kind of peculiarities observed in other 
spectral regimes and the quoted detection rates are to be considered 
as upper limits and are probably overestimated. 

As indicated in section 4.4, 487 of the 614 SCORPIO sources 
have a counterpart in at least one of the inspected catalogues. In 
particular, for our purposes, it is important to analyse how many 
SCORPIO sources have matches in stellar catalogues. The NO¬ 
MAD (Naval Observatory Merged Astrometric Dataset) catalogue 
is a major stellar catalog which contains data from the Hippar- 
cos, Tycho-2, UCAC-2 and USNO-Bl catalogues and is supple¬ 
mented by photometric information from the 2MASS final release 
point source catalogue. As indicated in Table IH most of the SCOR¬ 
PIO sources found in NOMAD do not have an optical counterpart 
(301/320 NOMAD matches are actually from 2MASS). 

This may indicate that a significant fraction of the SCOR¬ 
PIO sources represents a reddened stellar/galactic population. It 
is possible that this reddened population has some contaminants, 
namely some background galaxies. In principle, we may use the 
WISE colours to separate extragalactic versus Galactic sources, 
as th ose have revealed to be very effective in extragalactic sur¬ 
veys jjarrett et al.ll201lh . Stars are characterized by WISE colours 
clustering around zero magnitude while different classes of extra¬ 
galactic objects are clearly separated from them because of their 
redder colours. However, there are some Galactic objects, with 
extended dusty envelopes (AGB stars, YSOs and PNe) that have 
WISE colours similar to those of the different extragalactic popu¬ 
lations and they occupy th e same regions in the colour-colour dia¬ 
gram (iNikutta et al.ll20l4) . Therefore, it appears that, in our case, 
the WISE colour-colour diagram is of limited use for pointing out 
possible extragalactic contaminants and further ancillary observa¬ 
tions are necessary. 

Our observations, with unprecedented sensitivity and angular 
resolution at 2.1 GHz, have allowed, for the first time, to point out 
radio emission associated to some extended sources, many of which 
are classified as Galactic bubbles. Eor some of them the SCOR¬ 
PIO map is the first radio detection, while for the sources already 
known to be associated with a radio source, the SCORPIO map 
constitutes an impressive improvement with respect previous ob¬ 
servations. Eurthermore, comparison with existing mid-IR/FIR im¬ 
ages has helped to classify some of the objects. 
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Figure 9. Composite picture of the field around the position of IRAS 
16520-4146. The mid-IR/FIR maps from Spitzer (IRAC, 8 pm, green) and 
Herschel (PACS, 70 pm, red) are superimposed on the SCORPIO map 
(blue) with colour coding as indicated. In the insert in the low right cor¬ 
ner only the SCORPIO map (blue) is shown. 


REFERENCES 

Abbott D.C., Bieging J.H., Churchwell E., Torress A.V., 1986, 
ApJ, 303, 239 

AMI Consortium, Eranzen T. M. O., Davies M. L., et al., 2011, 
MNRAS, 415, 2699 

Benjamin R. A., Churchwell E., Babler B. L. et al., 2003, PASP, 
115,953. 

Beltran M. T., Brand J., Cesaroni R., et al., 2006, A&A, 447, 221 
Bieging J. H., et al., 1982, ApJ, 263, 207 
Bieging J. H., Abbott D. C. & Churchwell E. B. 1989, ApJ, 340, 
518 

Carey S. J., Noriega-Crespo A., Mizuno D. R., et al., 2009, PASP, 
121,76 

Churchwell E., Povich M. S., Allen D., et al., 2006, ApJ, 649, 759 
Chu Y.-H., Gruendl R. A., Guerrero M. A., et al., 2009, AJ, 138, 
691 

Condon J. J., Cotton W. D., Greisen E. W, et al., 1998, AJ, 115, 
1693 

Contreras Y, et al., 2013, A&A, 549, A45 
Cornwell T. J., 2008, IEEE Journal of Selected Topics in Signal 
Processing, Vol. 2, Issue 5, p.793-801 
Cutri, R.M. et al. 2014, VizieR Online Data Catalog, 2328, 0 
Cutri R.M. et al., 2003, VizieR Online Data Catalog, 2246, 0 
Drake S.A., Simon T., Linsky J.L., 1989, ApJS, 71, 905 
Dulk G. A. 1985, ARAA, 23, 169 

Dutra C. M., Bica E., Soares J., Barbuy, B., 2003, A&A, 400, 
533D 

Egan M.P., et al., 2003, VizieR Online Data Catalog, 5114, 0 
Fomalont E.B., 1989, in Perley R.A., Schwab F.R. & Bridle A.H., 
eds. Synthesis imaging in Radio Astronomy. A.S.P. CONE. SER. 
Vol 6, p. 213 

Fontani R, Beltran M. T., Brand J., et al., 2005, A&A 432, 921 
Eranzen T. M. O., Sadler E. M., Chhetri R., et al., 2011, MNRAS, 
439, 1212 

Gudel M. 2002, ARAA, 40, 217 

Gvaramadze V. V, Kniazev A. Y. & Fabrika, S. 2010, MNRAS, 
405, 1047 

Hales C. A., Norris R. R, Gaensler B. M., et al. 2014, MNRAS, 


© 2015 RAS, MNRAS 000, [T]-?? 















SCORPIO: A deep radio survey from the stellar life-cycle 11 


441, 2555 

Hallinan G., Antonova A., Doyle J.G., Bourke S., Lane C., Golden 
A., 2008, ApJ 684, 644 

Helfand D. J., Schnee S., Becker R. H., White R. L. & McMahon 
R. G., 1999, AJ, 117, 1568 

Helou G. & Walker D. W., eds. 1988, Infrared astronomical satel¬ 
lite (IRAS) catalogs and atlases. Volume 7: The small scale struc¬ 
ture catalog, Vol. 7 
Hogbom J. A., 1974, A&AS,15, 417 
Hogg 1989, AJ, 98, 282 

Ingallinera A., Trigilio C., Umana G., et aL, 2014, MNRAS, 437, 
3626 

Ishihara D., et ah, 2010, A&A, 514, A1 

Jarrett T. H., Cohen M., Masci K, et ah, 2011, ApJ , 735, 112 

Kimball A. E., Knapp G. R., Ivezic Z., et ah, 2009, ApJ, 701, 535 

Kronberger M., et aL, 2006, A&A, 447, 921 

Leitheree et aL, 1991, ApJ, 377, 629 

MacLeod G. C., van der Walt D. J., North A., et aL, 1998, AJ, 116, 
2936 

Mannings V, 1995, Ap&SS, 224, 185 
Mercer E. R, et aL, 2005, ApJ 635 560 

Mizuno, D. R., Kraemer, K. E., Flagey, N., et aL 2010, AJ, 139, 
1542 

Molinari S., Swinyard B., Bally J., et aL, 2010, PASP, 122, 314 
Murphy T, Mauch T., Green A., Hunstead R.W., Piestrzynska B., 
Kels A.P, Sztajer P 2007, MNRAS, 382, 382 
Nikutta R., Hunt-Walker N., Nenkova M., Ivezic Z. & Elitzur, M., 
2014, MNRAS, 442, 3361 

Norris R. P, Afonso J., Appleton P N., et aL, 2006, AJ, 132, 2409 
Norris R. P, Hopkins A. M., Afonso J., et aL, 2011, PASA, 28, 
215 

Osten R. A., Brown A., Ayres T. R., et aL, 2004, ApJS 153, 317 
Perry C. L., Hill G., Younger P. F. & Barnes J. V, 1990, A&AS, 
86, 415 

Purcell C. R., Hoare M. G., et aL, 2013, ApJSS, 205, 1 
Reynolds J., 1994, AT/39.3/040 

Rich J. W, de Blok W. J. G., Cornwell T. J., Brinks E., Walter K, 
Bagetakos L, Kennicutt R. C., Jr., 2008, AJ, 136, 2897 
Route M. & Wolszczan A., 2012, ApJL 747, L22 
Sana H., Rauw G., Naze Y, Gosset E., Vreux J.-M., 2006, MN¬ 
RAS, 372, 661 

Sanchez-Monge A., Beltran M. T., Cesaroni R., et aL, 2013, A&A, 
550, A21 

Seaquist E.R., 1993, RPPh, 56, 1145 

Setia Gunawan D. Y. A., Chapman J. M., Stevens 1. R., Rauw G., 
Leitherer C., 2003, lAUS, 212, 230 
Skrutskie M. K, Cutri R. M., Stiening R., et aL, 2006, AJ 131, 
1163 

Slee O.B., Wilson W, Ramsay G., 2008, PASA, 25, 94 
Soldi S., Shaw S., Miller J., et al. 2007, The Astronomer’s Tele¬ 
gram, 987, 1 

Spitzer Science Center, 2009, VizieR On-line Data Catalog, 11/293 
Sung H., Bessell M.S., Lee S.W 1998, AJ, 115, 734 
Trigilio C., Leto P, Umana G., Leone R, Buemi C.S., A&A, 2004, 
418, 593 

Trigilio C., Leto P, Umana G., Buemi C.S., Leone R, 2008, MN¬ 
RAS, 384, 1437 

Trigilio C., Leto P, Umana G., Buemi C.S., Leone R, 2011, MN¬ 
RAS, 739, LIO 

Umana G., Trigilio C., Catalano S., 1998, A&A, 329, 1010 
Umana G., Trigilio C., Hjellming R.M., Catalano S., Rodono M., 
1993, A&A, 267, 126 


Umana G., Trigilio C., Tumino M., Catalano S., Rodono M., 1995, 
A&A, 298, 143 

Urquhart J.S., et al., 2007, 461, 11 

Wachter S., Mauerhan J. C., Van Dyk S. D., et al., 2010, AJ, 139, 
2330 

Walsh A. J., Burton M. G., Hyland A. R. & Robinson G., 1998, 
MNRAS, 301, 640 

Walsh A. J., Hyland A. R., Robinson G. & Burton M. G., 1997, 
MNRAS, 291, 261 

Walsh A. J., Purcell C. R., Longmore S. N., et al., 2014, MNRAS, 
442, 2240 

White R. L., Becker R. H. & Helfand D. J., 2005, AJ 130, 586 
White S. M. 2000, Radio interferometry: the saga and the science, 
86 

White S. M., 2004, New Astronomy Reviews, 48, 1319 
White R.L., Becker R.H., Helfand D.J. 2005, AJ, 130, 586 
Wieringa M.H., Kesteven M.J, 1992, Measurements of the ATCA 
primary beam, AT Memo 39.3/024 
Wilson W. E., et al., 2011, MNRAS, 416, 832 
Wright A. E., Griffith M. R., Burke B. R & Ekers R. D., 1994, 
ApJS 91, 111 

Wood D. O. S. & Churchwell E., 1989, ApJS 69, 831 
Yamamura L, et al., 2010, VizieR Online Data Catalog, 2298, 0 
Zacharias N., Monet D.G., Levine S.E., Urban S.E., Gaume R., 
Wycoff G.L., 2004, AAS, 205, 4815 
Zoonematkermani S., Helfand D. J., Becker R. H., White R. L. & 
Perley R. A., 1990, ApJS, 74, 181 


© 2015 RAS, MNRAS 000, [I]-?? 


